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Abstract 

Single-molecule experiments provide new insights into biological processes 
hitherto not accessible by measurements performed on bulk systems. We report 
on a study of the kinetics of a triple-branch DNA molecule with four conforma- 
tional states by pulling experiments with optical tweezers and theoretical mod- 
elling. Three distinct force rips associated with different transitions between the 
conformational states are observed in the folding and unfolding trajectories. By 
applying transition rate theory to a free energy model of the molecule, probability 
distributions for the first rupture forces of the different transitions are calculated. 
Good agreement of the theoretical predictions with the experimental findings is 
achieved. Furthermore, due to our specific design of the molecule, we found a 
useful method to identify permanently frayed molecules by estimating the number 
of opened basepairs from the measured force jump values. 

Keywords: Nonequilibrium systems, single-molecule experiments, 
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1 Introduction 

In recent years, single-molecule experiments became of great importance in biophysi- 
cal research since progress in nano- and microscale manufacturing technologies facili- 
tated the design of scientific instruments with sufficient sensitivity and precision to en- 
able the controlled manipulation of individual molecules (for reviews, see, for example, 
l22l[T()lD . In contrast to the traditionally used bulk assays, where individual biomolec- 
ular dynamics can get masked, single-molecule experiments provide new insights into 
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the thermodynamics and kinetics of biophysical and biochemical processes hitherto not 
accessible. They complement standard spectroscopy and microscopy methods used in 
molecular biology and biochemistry and hence have to be regarded as an important 
source of additional information helping in the interpretation of biomolecular pro- 
cesses. Furthermore, single-molecule experiments permit the measurement of small 
energies and the detection of large fluctuations. 

The manipulation of single molecules offers a powerful new tool in molecular and 
cellular biophysics allowing for the exploration of processes occuring inside the cell at 
an unprecedented level. To instance just a few of the recently investigated biochem- 
ical processes: the transport of matter through pores or channels lfl9l l20l HI. interac- 
tions between DNA and proteins lfl3l or DNA and RNA ll36ll . the motion of single- 
molecular motors [2 1 4, 32], DNA transcription and replication l34lfT31 . virus infection 
ll26l 181. DNA condensation 11231 and ATP generation 11351 . In addition, the structure 
of biological networks l30l and the viscoelastic and rheological properties of the DNA 
(27] ED 0|29] have been studied. 

An important class of single-molecule experiments are performed with optical tweez- 
ers. By means of an optical trap generated by a focused laser beam, this useful tech- 
nique renders it possible to exert forces on micron sized objects, achieving sub-piconewton 
and sub-nanometer resolution in force and extension, respectively. Accordingly, one 
can study force-induced folding-unfolding dynamics and in this way get insight into 
corresponding processes in the cell and typical bond forces. Of particular interest is 
the unfolding of DNA molecules, where the hydrogen bonds between the complemen- 
tary base pairs (bps) are disrupted. This so-called unzipping is connected to the DNA 
replication mechanism. 

An interesting field of biophysical studies is the investigation of junctions in molecules 
since they present manifold ways to interact with other substances, for instance cations. 
Three-way junctions are especially interesting because metal ions such as magnesium 
can bind to them and alter the tertiary structure. Here a first step of such a study is 
presented where we investigate a molecule with a three-way junction alone, without 
cation binding. 

Many of the single-molecule experiments so far focused on molecules with a rel- 
atively simple free energy lanscape (FEL) exhibiting just two states, a folded and an 
unfolded one, or including an additional misfolding state, leading to different kinetic 
pathways. In this work we will consider a richer situation, where metastable states as 
intermediates occur during the folding-unfolding route. In this context, we will address 
the following key questions: 

(i) Can a corresponding molecule with such intermediate states be designed on the 
basis of a suitable model for a FEL? 

(ii) Is it possible to observe the intermediate states by perfoming pulling experiments 
with optical tweezers? 

(iii) Can phenomenological Bell-Evans kinetic models be applied to describe the 
folding-unfolding processes including intermediate states? In particular, when 
validating the kinetic theory against the experimental results, how do the first 
rupture force distributions compare with the ones predicted by the theory? 

A further important aspect that we looked at in some detail is the heterogenity 
of molecular folding-unfolding behaviour that we observed in the experiments. By 
measuring force-distance curves (FDCs) of several molecules we classified them into 
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different reproducible patterns. This leads to a useful method to identify irreversible 
molecular fraying, a phenomenon which is often observed in single-molecule studies. 



2 Description of the experiments 

Based on Mfold folding predictions ll25l [371 and taking FEL considerations into ac- 
count (see sec. [5]l, we designed and synthesised a DNA molecule which is composed 
of three parts and hence referred to as triple-branch molecule. It consists of a stem 
as introduced in (2JI with 21bps and two nearly identical hairpin branches which 
are formed of 16 bps and a loop with four bases, thus comprising a total number of 
114 bases, see fig. [T] To avoid misfolding, the second hairpin branch differs at two 
positions from the first one. 




Figure 1: Structure of the triple-branch DNA molecule. 

The triple-branch molecule is inserted between two identical short double-stranded 
DNA (dsDNA) handles of 29 bps each |9), leading to a total number of 172 bases, 
corresponding to a total contour length in the unfolded state of about lOOnm. Each of 
these polymer spacers is chemically linked to a bead. One of the beads is retained with 
the help of a pipette via air suction, the other one is optically trapped in a laser focus 
11281 . On the 5' end of the DNA molecule, biotin is attached to enable a connection with 
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a streptavidin-coated bead (SA bead) whose diameter is 1.8 \im. Biotin is a vitamin 
which establishes a strong linkage to the proteins avidin and streptavidin. The 3' end is 
modified with the antigen digoxigenin, able to interact with an antidigoxigenin-coated 
bead (AD bead). The latter has a diameter of 3.0 \xm. Figure [2] shows the different 
components of the molecular construct, that is to say the triple-branch DNA molecule, 
the handles and the beads, captured in optical trap and micropipette, respectively. Note 
that it is not a true-to-scale representation. 

The pulling experiments are carried out with a miniaturised dual-beam laser optical 
tweezers apparatus IfTTI at room temperature 25 °C) and at salt concentration of 1 
M NaCl aqueous buffer with neutral pH (7.5) stabilised by Tris HC1 and 1 M EDTA. 
The dual-beam optical tweezers collect data at 4 kHz and can operate with a feedback 
rate of 1 kHz. Spatial resolution constitutes 0.5 nm with a maximal distance range of 
^10 ^m. Forces up to 100 pN can be achieved, whereas the force resolution is 0.05 pN. 



DNA triple-branch molecule 




Figure 2: Sketch of the experimental setup. In our experiments we measure the relative 
distance X rather than its absolute value. 

Optical tweezer pulling experiments permit the measurement of the force / as well 
as the total distance X between the centre of the optical trap and the tip of the mi- 
cropipette, see fig. [2] In the experiment we vary the trap-pipette distance X(t) with 
a constant speed v = dX/dt in the range of 45 nm/s to 200nm/s, which corresponds 
to a constant average loading rate r of 3.0pN/s and 13.4pN/s in between rip events, 
respectively. The experiment consists of loading cycles which in turn are divided into 
an unfolding part (during "pulling") and a folding part (during "pushing"). The load- 
ing cycles are repeated as long as the tether connection is unbroken. Otherwise a new 
connection has to be established, possibly a new molecule must be searched and linked 
to a new bead. In sec.|4]we will discuss different patterns found in the measured curves 
and present a detailed analysis of two representative molecules. We chose them among 
seven molecules exhibiting the first and among five molecules featuring the second 
pattern. For each molecule we recorded, on average, approximately 50 cycles. The 
pulling speeds, ranging from 45 to 200 nm/s, influence the experimental results only 
weekly due to a logarithmic dependence of the first rupture force with the speed. We 
found compatible data for sets of similar molecules. In the theoretical analysis of the 
data in sec. [5] we concentrate on the largest set of 82 loading cycles for the molecule 
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at a speed of 200 nm/ s. The study of the above mentioned second molecule comprises 
55 cycles. 

3 Analysis of unfolding and folding trajectories 

Based on the design of the triple-branch molecule, we have to distinguish between four 
conformational states (see fig. [3]): 

1. a completely folded molecule. 

2. a completely unfolded stem with the hairpin branches still folded. 

3. stem and either hairpin branch 1 or 2 are completely unfolded. 

4. a completely unfolded molecule. 



1 . 2 . 3 . 4 




/ = 2nd + nioopd 



Figure 3: The four stable or metastable states of the triple-branch molecule: 1 - folded 
molecule, 2 - unfolded stem, 3 - stem and one hairpin branch unfolded, 4 - unfolded 
molecule. The values of x(n, f) refer to the end-to-end distance given in eq. with 
the number n of opened bps corresponding to the molecular construct shown in fig.Q] 
and the contour length / according to eq. (0. 

Figure [4(a)| displays a typical unfolding route and fig. |4(b)| a typical unfolding and 
folding trajectory in form of a FDC, which in fact records the evolution of the force as a 
function of time and relative trap position With rising X, the force / first increases 
almost linearly according to an elastic response of the DNA handles, which consist of 
dsDNA and are stable over the whole range of forces where the unfolding / folding 
of the triple-branch molecule takes place. The overstretching transition of the linkers, 
typically at 65 pN, lies much above the forces we explore. At a first rupture force f\ 
a sudden decrease ("jump") A/i occurs, which is caused by the unfolding of the stem. 
This unfolding goes along with an abrupt change in the length of the molecule when 
single-stranded DNA (ssDNA) is released. As a consequence, the bead in the optical 

1 In our experiments we measure the relative distance between trap and pipette, X, rather than the absolute 
value. The force / exerted on the molecular construct leads to a displacement f/k\, of the bead in the optical 
trap, where 1% = 0.08 pN/nm is the rigidity of the trap. Hence, the relative distance X is related to the relative 
molecular extension x m (see fig. [3) by x m = X — f/k^. 
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Figure 4: Force as a function of change in the trap-pipette distance (a) for one typical 
unfolding trajectory for the first investigated molecule and (b) for one typical unfolding 
and folding trajectory for the second investigated triple-branch molecule. Indicated are 
the three first rupture forces belonging to the transitions 1 (f\), 2 (fi) and 3 (/j) and the 
respective refolding forces, labelled fi^. The inset in (a) shows four further unfolding 
curves for the first investigated molecule. 

trap moves towards the centre of the trap, visible as the force drops. Following the 
jump A/i , there is again a linear increase up to the next force rip at a first rupture force 
fz, where one of the hairpin branches unfolds, which in turn leads to the force jump 
A/i- Eventually, the second hairpin branch unzips at a first rupture force fa with a jump 
A/3. The linear regime following this last force rip corresponds to the stretching of the 
whole molecular construct including handles and the already unfolded triple-branch 
molecule. 

Upon decreasing X from the completely unfolded state 4, the force first follows 
closely the corresponding unfolding part of the trajectory. However, a backward tran- 
sition does not occur at (/j - A/3) =; 16.2 pN, but at a considerably lower value of 
(14.0 + 0.8) pN, cp. fig. |4(b)| hence manifesting a hysteresis effect. Moreover, an in- 
vestigation of a larger number of folding trajectories reveals that the bases of the hairpin 
branches do not always pair conjointly in well-defined events during a short time inter- 
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val. In contrast, the corresponding force rips during unfolding indicate a cooperative 
behaviour of the biomolecule, where the breakage of all hydrogen bonds stabilising the 
DNA structure happens almost simultaneously. The second transition takes place at 
(13.3 ± 0.8) pN instead of (f 2 - A/ 2 ) — 15.5 pN. After both hairpin branches refolded, 
one can identify another sharp transition to the folded state 1 around (1 1.2 + 0.8) pN, 
which is again considerably lower than (f\ - A/i) ^ 13.7 pN. 

During unfolding, one can assume that the breakage of hydrogen bonds follows the 
sequence of the molecular construct. This makes it useful to introduce the number n 
of broken bonds as state variable and to calculate a FEL as function of this variable 
(see sec. 0. The refolding of the hairpin branches in the folding trajectories exhibit 
less sharp transitions, see fig. |4(b)| During folding, in particular at the beginning in the 
unfolded state, a huge number of secondary structures can be found which implies that 
the kinetic pathways are less predefined and accordingly, the transitions get smeared 
out. With respect to a theoretical treatment, moreover, a description in terms of the 
simple state variable n becomes unlikely to be sufficient. In a refined analysis, many 
more configurations should have to be included as relevant states in a coarse-grained 
description [171 . Such refined analysis, however, goes beyond the scope of this work 
amd we therefore concentrate on the unfolding process in the following. 

There are plenty of possible ways to analyse the FDCs in order to find out the first 
rupture forces and the force jumps. In our procedure we arranged the normalised data, 
i.e. the relative distance X and the force /, in windows of a certain size of data points. 
For all consecutive windows we then calculated the slope of the considered data points, 
the span, i.e. the maximum distance between the lowest and the highest force value, 
and the mean of the relative distance Xj as a moving average. Transitions between 
the conformational states take place where the slope is minimal and the span maximal 
under the condition that an appropriate number of contiguous windows is connected. 
Having found the Xj of the three force rips, slope and axis intercept are calculated by 
linear regression for each conformational state. One can now easily calculate the first 
rupture forces as the intersection points with the four fitted lines and extract the force 
jump values, as exemplified in fig. [5] for both molecules whose unfolding trajectories 
were depicted in figs. |4(a)| and |4(b)| respectively. Since the data acquisition rate is 




Figure 5: The first rupture forces and force jump values of the unfolding trajectories 
shown in fig. [4] are extracted as indicated here from the intersection of the fitted grey 
lines and the transitions (vertical lines). In part (a) the procedure is shown for the 
molecule used in the theoretical analysis in sec. [5] and in (b) for the molecule with 
permanent fraying behaviour. 
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constant, at higher pulling speeds less data points are collected and therefore the values 
of the analysis are broader distributed. The three first rupture forces fi as well as 
the jumps Afi are subject to stochastic fluctuations, as can be seen in fig. |4(a)| where 
we show four unfolding trajectories belonging to different pulling cycles of the same 
molecule in the inset. An analysis revealed that the fluctuations of the force jumps Afi 
are about ten times smaller than the fluctuations of the /;. Accordingly, in sec. [5] we 
will disregard the fluctuations in the A/; and use only their averages Af, that will be 
discussed in more detail in the following section. 

4 Different unfolding patterns related to the number of 
opened bps 

Applying the above mentioned procedure to analyse the experimental data of several 
molecules, we found two predominating patterns in the unfolding trajectories which 
are reflected in the distributions of the first rupture forces f, as follows. In the first 
pattern, see fig. |6(a)| these distributions have a similar shape for all three force rips. The 
histograms indicate the existence of one maximum slightly below 17 pN. The mean 
values for the three first rupture forces are f\ = (17.0 + 0.9) pN, f 2 = (16.4 ± 0.6) pN 
and fi = (16.8 + 0.7)pN. In contrast, in a second pattern we detected a strikingly 
lower value for the first rupture force of the first rip f\ — (14.6 + 0.8) pN, as depicted in 
fig. |6(b)| whereas the other two rip forces lie basically in the same range of 16 to 18 pN. 
As in the former case, the second rip tends to have a slightly smaller first rupture force, 
J 2 = (16.5 + 0.8) pN, than the third rip, % = (17.2 + 0.5) pN. 

The small f\ observed in the second molecule suggests the occurrence of permanent 
molecular fraying. Obviously less force is needed to unfold the stem than typically, cp. 
fig. |4(a)| with |4(b)| since some bps at its basis are partly or completely melted. In other 
words, during folding, this molecule does not reach an entirely folded state but some 
bps of the stem next to the handles remain irreversibly and permanently open. This 
phenomenon has been observed previously in several pulling experiments lfTTl[33ll . A 
possible reason for irreversible fraying is the formation of reactive oxidative species 
due to the impact of the laser light of the optical trap leading to a degradation of the 
DNA bases 1121 . The so generated singlet oxygens are known to oxidise certain nucleic 
acids, such as guanine and thymine, irreversibly. This could explain our observation 
that once a molecule shows fraying it does not change back again to normal behaviour. 
Due to the fact that we work with polystyrene microspheres which are more prone to 
photodamage than the DNA bases themselves, their wide ranging interaction with the 
bases might be reduced replacing polystyrene by silica beads, which exhibit consid- 
erably minor irreversible oxidative damage. It would be very interesting to carry out 
such experiments. 

The insets in figs. |6(a)| and |6(b)| depict the corresponding force jump distributions 
of both molecules. While the first molecule possesses a large first force jump value 
of A/i = (1.3 + 0.2) pN and two smaller force jumps at the second and third rip of 
Af 2 = (0.9 ± 0.2) pN and A/ 3 = (1.0±0.2)pN, the frayed molecule features three 
force jumps of approximately the same value, i.e. A/i = (0.93 + 0.05) pN, Af 2 = 
(0.95 + 0.06) pN and A/ 3 = (0.97 + 0.04) pN, respectively. This illustrates clearly the 
influence of irreversible fraying in the latter case since, roughly estimated, the same 
number of bps is expected to open in all three rips which should not be the case in an 
entirely folded molecule. 
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Figure 6: Histograms of the three first rupture forces during unfolding for two repre- 
sentative triple-branch molecules: (a) the one used in the theoretical analysis in sec. [5] 
and (b) the one exhibiting permanent fraying behaviour. Note that the average force 
value of the first rip has decreased in (b) as compared to (a). The insets show the 
corresponding histograms for the force jump values. 

During a force rip, the relative distance X is constant and thus AX = 0. Therefore 
the change in the relative molecular extension, 

Ax m = AX - Af/kes = Ax(n, /), ( 1 ) 

is only related to the force jump Af and the combined stiffness of bead and handles 
^df, given by l/k^ = l/k\, + l/k^, where k\, is the trap stiffness and the rigidity 
of the handles, respectively. With the help of a linear least squares fit, the effective 
stiffness k e s of the molecular construct is extracted from the average slope of the FDCs 
and amounts to (0.067 + 0.003) pN/nm. Note that, as expected, k e $ is smaller than k\, 
O0.08pN/nm). 

For a certain force value / and assuming an elastic model for the released ssDNA, 
the number n of opened bps is related univocally to the equilibrium end-to-end distance 
of the DNA molecule x(n,f), whereas its change Ax(n,f), in turn, equals Ax m . Using 
this relation, it is now possible to estimate the change in the number An, of bps which 
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are opened sequentially during each force rip. Considering only the configurations of 
the four conformational states, the equilibrium end-to-end distance x(n,f) can be de- 
composed into two parts, cp. fig. [3] The first part, the elongation ui(f) of the mean 
end-to-end distance of the ssDNA along the force direction, accounts for the ideal elas- 
tic response of the ssDNA, where / is the contour length. The second part contains the 
contribution of the diameter of stem and hairpin branches, respectively. Acccordingly, 

( 0, n =21 + 16+ 16 bps 
x(n,f) = u,(f) + I d , n = and 21 + 16 bps , (2) 
[ d' , n =21 bps 

where do ^ 2 nm, in accordance with the diameter of the B-DNA helix. The exact value 
of the diameter contribution of both hairpin branches dL when the stem is unfolded, 
depends on the orientation of the branches. We set d' Q - 2do as a working valu^. 

Regarding the contour length /, which depends, amongst others (cp. sec. [5), on the 
number n of opened bps, and considering again solely the configurations of the four 
conformational states, one gets 

( 0, n = 0and21bps 

/ = 2nd + < n\ oop d, n =21 + 16 bps , (3) 
{ 2ni mp d, n =21 + 16 + 16 bps 

where the interphosphate distance d is taken to be 0.59 nm/base and «i oop = 4 is the 
number of bases per loop. For every force rip f, the corresponding number of opened 
bps «, is calculated separately by considering the differences of / between the states. 

Different types of models can be used to calculate m;(/). Prominent examples bor- 
rowed from polymer physics are the freely jointed chain (FJC) and the worm-like chain 
(WLC) model. According to ref. 1271 . the FJC model includes an extra term, leading 
to the expression 

bf\ k B T 



coth 



k B Tj bf 



(4) 



Here Y denotes the Young modulus, b is the Kuhn length, k B the Boltzmann constant 
and T the temperature. Typical values of the Kuhn length and the Young modulus under 
working conditions of T - 25 °C and 1 M NaCl concentration are b = 1 .42 nm and 
Y = 812pN 11271 or, as published recently, b = 1.15 nm and Y = °o IfTTI . respectively. 
In the WLC model (3), the force f(ui), due to an elongation ui, is given by 



k B T 



1 1 Ui 



4(1 - M/ //) 2 4 I 



(5) 



and to obtain «/(/), this equation has to be inverted. Based on the WLC model, we 
tested the influence of the persistence length P in a typical range of 1.0 to 1.5 nm f6l . 

From the data shown in table Q] it is apparent that, depending on the model and 
parameters, the results for the estimated change in the number of opened bps vary in an 
acceptable range when the errors are taken into account. In addition, one can see that it 
is not evident which model and parameters should be considered as the best ones. Good 
results are found for the FJC model using recent values of IfTTI and for the WLC model 



2 Due to this simplification, the A«i of the first rip is likely to be slightly underestimated and the second 
rip's A?i2 overestimated. However, it will not affect the change in the total number An tot of opened bps since 
we consider the change of x(n, /), and the d' contributions will cancel each other out. 
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molecule 1 


molecule 2 


model 


parameter 
b & P [nm], Y [pN] 


Awi 


An 2 


change 
Aw 3 


in no. of 

Antot 


opened b] 
Awi 


)s [bps] 

An 2 


An 3 


An tot 


FJC 


b = 1.42, Y = 812 |27| 
b = 1.15, 7 = oo rm 


18(2) 
19(3) 


14 (2) 
15(2) 


15(2) 
16(2) 


46 (3) 
50 (3) 


13(1) 
14(1) 


15(1) 
16(1) 


15(1) 
16(1) 


42 (2) 
46 (2) 


WLC 


P = 1.0 
P= 1.3 11241 
P= 1.5 


20 (3) 
20 (3) 
19(3) 


16(3) 
15(2) 
15(2) 


17(3) 
16(2) 
16(2) 


54 (3) 
51 (3) 
50 (3) 


14(1) 
14(1) 
13(1) 


17(1) 
16(1) 
16(1) 


17(1) 
16(1) 
16(1) 


49 (2) 
46(2) 
45 (2) 


expected values 


21 


16 


16 


53 


21 


16 


16 


53 



Table 1 : Overview over the change in the number of opened bps for different models 
and parameters for the two representative molecules. The numbers in brackets are the 
standard deviations. 



with P = 1.3 nm [24]. We chose to work with the FJC model with the parameters of 
IfTD in sec.El 

The second molecule indeed reveals a considerably smaller An i , depending on the 
model around 13 or 14 bps, so that 7 or 8 bps are not closed after the folding process 
is completed. Performing single-molecule experiments without knowing the exact in- 
fluence of permanently frayed bps can lead to misinterpreted results. Checking the 
appropriate parameters for the polymer models with the help of the change in the num- 
ber of opened bps of the hairpin branches, one can estimate the number of irreversibly 
frayed bps at the basis of the sternal 

5 Theory for the unfolding kinetics 

The kinetics of the unfolding process can be described on a coarse-grained level based 
on a Gibbs free energy G(n,f) as a function of the number n of sequentially opened 
bps for an applied force /. For small forces, including / = 0, the FEL is expected 
to have a shape as displayed in fig. |7(a)| In general, G(n, 0) increases monotonously 
with «. However, local minima occur at the metastable states 2, 3 and 4 because there 
is an increase of entropy associated with the release of additional degrees of freedom 
when the stem-hairpin-j unction and end-loops of the hairpin branches are opened. With 
rising force the FEL is expected to get tilted, so that the energies of the metastable states 
are lowered. With a knowledge of G(n, /) we can apply standard transition rate theory 
and write for the transition rate from state i to i + 1 

r w+ i(/) = 7- Yu+iCf), ( 6 ) 

where y° is an attempt rate and Yy+i(/) is the Boltzmann factor corresponding to the 
activation barrier AG,;, + i(/) that has to be surmounted, 

Yi,i+ 1 (/) = ex Pl -j—f I ■ ( 7 ) 

When considering only sequential configurations in the evaluation of the FEL, dif- 
ferent structures compatible with a given n can occur once the stem is completely un- 
folded. These refer to different possibilities of breaking the bps in the hairpin branches 



3 We like to note that the checking of the change in the number of opened bps can be, in principle, also 
applied to non-permanent, reversible molecular fraying. 
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Figure 7: (a) Sketch of the expected FEL for the triple-branch molecule at zero force 
as a function of the number « of opened bps. (b) FEL calculated from eqs. ([8]) to (flTT i 
at two different forces in units of the thermal energy k%T. 

1 and 2. In order to point out this "degeneration", a new parameter a is introduced, 
leading to the FEL G(n, a, /). At a given n, we must calculate the (restricted) partition 
sum over the configurations a to get G{n,f). In order to find G(«, a,/), we consider 
the following decomposition, 

G(«, a, f) = G foral («, a) + Gl(n, a, f) - fAxf , (8) 

where Gf 0lm (n, a) is the free energy of formation of the configuration (n, a), G" r (n, a, f) 
is the strain energy of the unfolded ssDNA and fAxf is a Legendre term (for the defi- 
nition of Axf see eq. ( TTOb below). 

The free energy of formation Gf orm («, a) is written as 

Gform(«, a) = ^ g^+i + G junc (n) + G, ( ^ p (n) + G® p (n) . (9) 

all bps 

The first term refers to the nearest neighbour model developed in |5]|7], which specifies 
the interaction g^+i between a base pair ju and the directly adjacent one fi + 1. It was 
shown to provide reasonable agreement with experiments lfl6l I2T1 l33l . For example, 
applying this model onto a sequence 5'-TCCAG. . . -3' and its complementary part 3'- 
AGGTC. . . -5', the stack energy reads G stac k = £tc/ag + gcc/GG + £ca/gt + £ag/tc + ■ • ■ 
The most recent values of g^+i lie in the range of -2.37 to -0.84kcal/mol at 25°C 
IfTTI . The terms Gj unc (n), G^ p (/i) and G^ p (/i) in eq. (0 describe the free energy re- 
duction due to the release of the stem-hairpin-j unction and end-loops and are estimated 
from |25] EH as G, ( ^ p = G® p = 1.58kcal/moland G jlmc = 4.90kcal/mol. 

The strain energy Gf tt (n, a, f) of the unfolded single-stranded part OTI with contour 
length I = /(«,o;jj can be calculated from the work needed to stretch the unpaired 
bases. We like to remind the reader that we denote the elongation of the mean end- 
to-end distance of the ssDNA in force direction by m/(/). In what follows we chose 

4 In previous publications of some of the authors this contour length was denoted by l„ a to emphasize the 
dependence on n (and, in addition, a here). This dependence is caused by the change of the contour length 
in the transitions. For easier reading we suppress to give it explicitely in the following. Further details about 
the contour length were already discussed in sec. [4] 
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to work with the FJC model, see eq. @, with parameters of 0T | instead of using the 
WLC model, despite the fact that both approaches give similar good results (cp. table 
[TJ. Since m/(/) is monotonously increasing with /, it has an inverse fi(u) = uj l (f), 
which is the force that is exerted by a ssDNA chain with contour length /, if its mean 
end-to-end distance is elongated by u. Accordingly, setting Axf = «/(/), we can write 

Axf / 

Gl(n,a,f)= J du'Mu') = f^- J df'uAf). (10) 
o o 

Finally, we computed G(n, f) by 

G{n,f) = -fc B rin^exp|- ^J^ j . (11) 

In fig. [7(b)] the FEL is depicted for / = and / = 16.21 pN. It exhibits the 
behaviour anticipated in fig. |7(a)| for zero force it has minima at the stable/metastable 
states and it becomes tilted with rising force. When approaching the force regime 
where the rips occur in fig. |4(a)| the levels of the minima become comparable. We want 
to point out that the DNA sequences shown in fig.[T]have been designed deliberately to 
yield the multiple-state structure seen in fig. |7(b)| This gives us some confidence in the 
model underlying the construction of the G(n, a, f) in eq. ([8]). 

Based on the FEL we can easily calculate the transition probability for 
the first rupture force f; in the ;th transition if the first rupture force was /;_] in the 
(z - l)th transition. The result is 

Wil/i-i) = - Ym+i(/) exp [-£ [ f df Ym + i(/')1 , (12) 

where f* = fi — A/j (and /o = / * = 0); r was the loading rate, see sec. [2] and y° and 
Yu+i(/) were defined in eq. ©. The activation energy AG i(i+ i(/) appearing in eq. (O 
was calculated, as indicated in fig. |7(a)| from the G{n,f) by determining the energy 
G™ n (/) of the local minimum belonging to state i and the saddle point energy G^'J < j le (/) 
of the rth transition between the ith and (/+ l)th state, AG i>M (f) = G^f e (/)-G™ n (/). 
The attempt rate y° was used as the only fitting parameter. 

For the joint probability density of the three first rupture forces we then obtain 

Wi./z./a) = WilO) W(m) W(/ 3 |/ 2 ) , (13) 

which allows us to calculate the distributions shown in fig. |6(a)| 

Figure |8] displays the histograms for the three rips from fig. |6(a)| in comparison 
with the distributions calculated from our theory. In view of the available statistics (82 
cycles, see sec.|2]i, the agreement is quite satisfactory. 



6 Conclusions 

An important class of biophysical studies is the investigation of junctions in molecules 
since they present manifold ways to interact with other substances, for instance cations. 
Three-way junctions are especially interesting because metal ions such as magnesium 
Ifl4l [181 can bind to them and alter the tertiary structure. Here a first step of such a 
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Figure 8: Comparison of the first rupture force distributions from fig. |6(a)| with the 
theory for the first (/i), the second (/i) and the third (fi) transition, see left, middle and 
right panel, respectively. The best fit was obtained for attempt rates y? = 9 • 10 6 Hz, 
y° = l.l-10 6 HzandY? = 3.5-10 5 Hz. 



study is presented where we investigate a molecule with a three-way junction alone, 
without in vivo relevant substances. 

One of our aims of this work was to study whether the construction and kinetics 
of more complex DNA molecules with richer folding-unfolding behaviour can be de- 
scribed by proper extensions of theories developed successfully for two-state systems 
so far. Our results show that this is indeed possible, at least for the unfolding trajec- 
tories. A triple-branch molecule has been specifically designed to produce a four-state 
system based on a model for the free energy landscape. This design was successful 
and we were able to prove the existence of these states by the emergence of associated 
force rips in pulling experiments. The first rupture forces have been systematically 
recorded in these pulling experiments and their distributions have been calculated. A 
transition rate theory based on the free energy landscape was successful in describing 
these distributions. 

Two patterns have been found in the measured unfolding trajectories, one indicat- 
ing the anticipated unfolding behaviour and the other one pointing to the occurrence 
of irreversible molecular fraying. This characterisation was possible by connecting the 
extracted force jump values to the change in the number of opened bps at each transi- 
tion. For this estimation we tested the validity of two polymer models for the elastic 
response of ssDNA (FJC and WLC) and different sets of parameters in order to find the 
best agreement with the expected values. This analysis is useful to compare the elastic 
properties measured in DNA unzipping experiments with those obtained by stretching 
ssDNA polymers IfTTI . 

One class of molecules required a smaller force than anticipated to unfold the stem 
since some bps at its basis are partly or completely melted due to photodamaging. 
Permanent molecular fraying is an usually undesired, but frequent effect in single- 
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molecule studies and deserves special attention in order to reduce its distorting influ- 
ence on experimental results. It is important to find means to avoid irreversible fraying 
since not fully closed molecules change the measured unfolding-folding trajectories so 
that an average over all molecules, including permanently frayed ones, can lead to de- 
viations of the real values and to misinterpreted results. To improve the statistics of the 
results, it is therefore necessary to identify irreversibly frayed molecules and remove 
them from the analysis. Within our analysis, we found a useful method to identify 
permanent molecular fraying. With the appropriate parameters for the polymer models 
one can estimate the number of irreversibly frayed bps at the basis of the stem. Further 
experiments performed with silica beads instead of polystyrene microspheres could 
clarify under which conditions permanent molecular fraying can be decreased. 

All this knowledge paves the way to further interesting studies such as the refold- 
ing problem IfTTl or two topics which require further experimental research. Firstly, 
the binding of metal ions to the three-way junctions could be examined in order to 
find out about structural changes due to the formation of tertiary contacts and there- 
fore altered kinetics of the unfolding process fTSl . Secondly, the translocation motion 
of helicases that unwind dsDNA could be addressed as another interesting subject, in- 
cluding the investigation of how they move along bifurcation points. Eventually, the 
kinetic approach for the prediction of the probability distributions of the first rupture 
forces could be a well suited starting point for future in-depth modelling in this domain 
on a more microscopic basis. 
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